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ABSTRACT 


The  eflFects  of  ocean  surface  waves  on  optical  signal  fluctuations  and  optical  refi-action 
are  evaluated.  The  effects  on  signal  fluctuations  were  inferred  fi’om  the  ocean  surface  wave 
spectral  density  and  temperature  fluctuation  spectral  density.  The  effects  of  ocean  waves  on 
optical  refraction  are  manifested  in  the  near-surface  gradient  of  the  temperature  profile.  A 
modified  temperature  profile  incorporating  a  statistical  mixing  length  scale  is  proposed.  The 
research  points  to  the  importance  of  certain  parameters  that  should  be  considered  when 
optical  signals  skim  the  ocean  surface. 
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1.  INTRODUCTION 


A.  BACKGROUND 

The  importance  of  environmental  effects  on  communication  links  and  radar  systems 
has  been  recognized  for  the  last  few  decades.  Low  altitude  electromagnetic  wave  propagation 
over  the  ocean  is  strongly  dependent  on  atmospheric  conditions  and  sea  state.  It  is  essential 
to  accurately  characterize  atmospheric  turbulence  in  order  to  predict  the  propagation  of 
electromagnetic  signals.  For  a  signal  skimming  the  ocean  surface,  an  understanding  of  ocean 
surface  wave  behavior  should  be  given  equal  emphasis  in  electromagnetic  wave  propagation 
modeling.  Specifically,  it  is  necessary  to  relate  the  variation  of  the  index  of  refi-action  to 
atmospheric  turbulence  and  sea  state. 

The  refractive  index  at  optical  fi'equencies  can  be  approximated  by 


«-l  =  77.6(1 +7.52xl0■3A•2)(P/7)10■^ 


(1) 


where  P  is  the  atmospheric  pressure  in  milhbars,  T  is  the  temperature  in  Kelvin  and  A  is  the 
wavelength  of  light  in  pm.  To  determine  the  variation  in  the  index  of  refi-action,  it  is  sufficient 
to  know  the  variation  of  temperature  in  the  atmospheric  surface  layer.  In  the  surface  layer, 
the  vertical  profiles  of  wind  speed,  w,  temperature,  T  and  humidity,  q,  all  have  logarithmic 
forms  ^  =(w,  T,  q),  given  by  [Ref  1] 


1 


^  [ln(z/2o5)  -  il;^ (C)]/k  , 


(2) 


where  is  the  flux  scaling  parameter,  is  the  value  at  the  surface,  is  the  roughness 
length,  K  is  the  von  Karman  constant,  C  =  L  the  Monin  Obukov  length  and  the 
empirically  derived  stability  correction  functions. 

Equation  (2)  shows  that  the  model  results  in  large  refractive  index  gradients  in  regions 
close  to  the  sea  surface.  An  ej^eriment  conducted  by  the  Naval  Research  Laboratory  in  June 
1994  [Ref  2]  showed  that  the  result  of  the  large  index  gradient  caused  rays  from  the 
transmitter  to  be  bent  more  than  that  required  to  establish  optical  continuity  between  the 
transmitter  and  the  receiver. 

Unlike  the  atmospheric  boundary  layer  above  land,  the  boundary  layer  above  the  sea 
is  continuously  interacting  with  the  ocean  surface  waves.  The  temperature  profile  of  Equation 
(2)  derived  for  an  immovable  surface  may  not  be  valid  in  the  region  of  a  few  meters  above  the 
sea  surface.  Ejqjerimental  data  on  signal  heights  that  commensurate  with  ocean  wave  heights 
demonstrate  deviations  from  the  logarithmic  distribution  given  by  Equation  (2)  [Ref  3]. 

B.  OBJECTIVES  OF  THE  RESEARCH 

This  research  evaluates  the  effects  of  ocean  surface  waves  on  optical  fluctuations  and 
relates  the  optical  refraction  to  the  near-surface  gradient  of  temperature  and  ocean  surface 
waves.  Measurements  from  an  optical  link  are  used  to  infer  the  temperature  profile  and  hence 
the  refractive  index  profile.  These  are  used  to  assess  the  effect  of  the  ocean  waves.  Based  on 
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the  experimental  data,  a  proposed  temperature  profile  is  used  in  conjimction  with  a  ray  tracing 
program  to  validate  optical  continuity. 

C.  ANALYSIS  APPROACH 

The  research  will  first  examine  the  small-scale  atmospheric  fluctuations  due  to  the 
presence  of  the  ocean  surface  waves  on  the  optical  signal  fluctuations.  Computed  ocean 
surface  spectral  density  and  air  temperature  fluctuation  spectral  density  are  used  to  establish 
the  influence  of  the  ocean  surface  waves.  Having  established  the  influence  of  the  ocean  waves, 
modification  to  the  near-surface  mean  temperature  gradients  is  proposed.  This  is  done  by 
incorporating  a  random  scaling  length  to  account  for  the  random  motion  of  the  ocean  waves. 
With  the  proposed  temperature  profile,  a  ray  tracing  code  was  used  to  establish  optical 
continuity  between  the  transmitter  and  the  receiver. 

This  report  comprises  two  volumes.  Volume  I  contains  seven  chapters  covering  the 
following  details:  The  experimental  set-up  and  instrumented  flux  buoy  are  described  in 
Chapter  II.  Chapter  El  presents  the  ray  tracing  code  used  in  conjunction  with  the  proposed 
temperature  profile  to  vahdate  optical  continuity  between  the  transmitter  and  the  receiver. 
The  analysis  of  ocean  surface  wave  spectral  density  is  presented  in  Chapter  IV  and  the  bulk 
aerodynamic  method  used  to  analj^e  the  meteorological  data  is  covered  in  Chapter  V.  The 
results  of  data  analysis  are  presented  in  Chapter  VI.  Chapter  VII  summarizes  the  findings  and 
recommendations  of  this  research.  The  complete  set  of  data  plots  is  contained  in  Volume  H. 
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n.  EXPERIMENTAL  SET-UP 


A.  OPTICAL  LINK 

In  the  experiment,  an  optical  link  at  0.83  pm  was  set  up  from  Point  Pinos  at 
Lat/Long.  (36®  38’  10.93 1’N,  121®  55'  45.444"  W)  to  Marina  State  Beach  at  Lat./Long.  (36® 
42’  01.365"  N,  121®  48’  33.651"  W)  .  The  propagation  path  was  approximately  12.8  km  and 
the  bearing  of  the  transmitter  from  the  receiver  was  237®.  A  sketch  of  the  optical  path  is 
shown  in  Figure  1. 

The  transmitter  consists  of  a  0.83  pm  GaAs  laser  source  with  an  average  power  of 
3  5mW  and  intensity  modulated  at  4.2  kHz.  It  has  a  beam  aperture  of  50  mm  and  a  beam 
divergence  half  angle  of  1.3  mrad.  The  receiver  has  a  50  mm  aperture,  a  laser  line  filter  with 
a  bandwidth  of  40  nm  centered  at  0.83  pm  followed  by  an  Oriel  Instrument  Inc.  battery 
operated  ^con  photodetector  ^d  pre-amplifier  assembly.  Block  diagrams  of  the  transmitter 
and  receiver  are  shown  in  Figures  2  and  3  respectively.  During  the  experiment,  the  heights 
of  the  transmitter  and  receiver  were  varied  and  the  received  signal  recorded  on  tape. 

The  recorded  signal  was  demodulated  with  a  Stanford  Research  Systems  lock-in 
amplifier  with  sensitivity  set  at  50  mV,  time  constant  of  1  ms  and  internal  oscillator  set  at  4.2 
kHz.  The  output  of  the  lock-in  amplifier  was  connected  to  a  Labview  National  Instruments 
12  bit  analog  to  digital  card  plugged  into  a  Macintosh  computer.  The  analog  signal  was 
sampled  at  400  Hz.  A  total  of  41  signal  samples  covering  the  period  from  30  May  96  to  5  Jim 
96  were  taken. 


5 


Figure  1.  Optical  Path  across  the  Monterey  Bay 
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GaAs  Laser  Diode 


Figure  2.  Transmitter  Block  Diagram 
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Figure  3.  Receiver  Block  Diagram 


B.  FLUX  BUOY 


The  flux  buoy  was  deployed  at  approximately  the  mid  point  of  the  optical  link,  at 
Lat./Long.  (36°  40' N,  121°  52’  W)  where  the  ocean  depth  is  about  69.5  m.  The  flux  buoy  was 
deployed  to  record  meteorological  and  wave  data.  The  flux  buoy  system  contains  two 
separate  sensor  suites,  viz.,  a  mean  meteorology  package  and  a  high  frequency  sampling  ‘flux’ 
package.  Two  sets  of  data  were  collected;  the  high  frequency  data  (flux)  of  air  temperature 
and  wind  velocity  fluctuations  sampled  at  21  Hz,  wave  motion  data  sampled  at  10.5  Hz  and 
low  frequency  data  (mean  meteorology)  of  wind  speed,  wind  directions,  humidity,  pressure, 
air  and  sea  temperature  sampled  at  one  minute  interval.  A  picture  of  the  buoy  is  shown  in 
Figure  4. 

The  high  frequency  air  temperature  and  velocity  were  measured  by  the  Solent 
ultrasonic  anemometer  at  5.09m  above  mean  tide  level.  The  measurement  technique  is  based 
on  sound  pulses  transmitted  back  and  forth  between  transducers.  The  time  of  flight  in  each 
direction,  say  tj  and  t2  are  measured.  If  c  is  the  speed  of  sound,  x  the  distance  between  the 
transducers,  and  there  is  an  air  flow  v  along  the  line  of  the  transducers,  the  following 
relationships  are  readily  derived: 


h  - 


X 

(c+v)  ’ 


(4) 
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V  =  0.5  X 


(5) 


1  1 


and 


c  =  0.5  X 


1+1 

h  ^2 


(6) 


By  arranging  three  pairs  of  transducers  in  different  orientations,  the  direction  and  magnitude 
of  the  incident  air  flow  may  be  imambiguously  derived.  The  air  temperature  is  calculated  from 
the  speed  of  sound  as  follows  [Ref  4]: 


=  yRT(l+0.51q)  , 


(7) 


where  yR  is  403  m^  s'^  K*  for  air,  T  is  the  air  temperature  in  degrees  Kelvin  and  q  is  the 
specific  humidity  in  kg  of  water  vapor  per  kg  of  air. 

The  ocean  surface  wave  data  was  measured  by  the  Systron  Dormer  Motion  Pak.  The 
Motion  Pak  consists  of  three  accelerometers  measuring  the  accelerations  of  the  x,  y  and  z 
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directions  (sway,  surge  and  heave)  and  three  rate  gyros  measuring  the  angular  rates,  0,  (f)  and 
ij;  (pitch,  roll  and  yaw).  The  data  was  sampled  at  10.5  Hz.  The  Motion  Pak  uses  a  right- 
handed  coordinate  system  for  the  angular  rates  about  the  x,  y,  and  z  axis.  It  was  mounted  so 
that  the  x-axis  points  towards  the  port  side  of  the  buoy,  the  y-axis  points  towards  the  bow  and 
the  z-axis  downward.  The  z-axis  voltage  was  adjusted  so  that  when  the  buoy  is  at  rest  on  a 
flat  surface,  the  z-axis  linear  acceleration  voltage  output  read  zero  voltage. 

Three  wave  staffs  were  deployed  to  measure  the  high  frequency  ocean  surface  waves. 
The  three  wave  staffs  were  arranged  such  that  each  staff  was  located  at  the  vertex  of  an 
equilateral  triangle.  The  arrangement  allows  for  the  computation  of  the  wave  directional 
spectra.  Unfortunately,  wave  staff  number  2  was  not  functioning  properly  during  the 
experiment. 

The  buoy’s  heading  relative  to  magnetic  North  can  be  determined  from  the  three 
Syntron  compass  outputs.  The  magnetic  deviation  from  the  true  North  at  the  buoy  was  15® 
15'  East.  The  buoy’s  heading  relative  to  the  true  North  can  be  determined  by  subtracting  the 
deviation  from  the  buoy  magnetic  heading. 

The  low  frequency  wind  speed  was  measured  by  the  R.M.  Young  Wind  Monitor 
located  at  3.80  m  above  mean  tide  level.  The  wind  speed  sensor  is  a  helicoid  shape  propeller 
molded  of  polypropylene  plastic.  Threshold  sensitivity  of  the  propeller  is  1.0  m/s.  Rotation 
of  the  propeller  produces  a  linear  DC  voltage  signal  directly  proportional  to  the  wind  speed. 

The  Everest  IR  sensor  was  used  to  measure  the  sea  surface  temperature.  The  IR 
radiation  from  the  sea  surface  is  collected  by  an  optical  lens  system  and  focussed  on  an  IR 
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detector.  The  detector  in  turn  converts  the  IR  signal  to  an  electrical  signal  proportional  to  the 
sea  surface  temperature. 

The  low  frequency  temperature  and  humidity  were  measured  by  the  Rhotronics 
relative  humidity  and  temperature  sensor.  The  humidity  sensor  is  a  small  hygroscopic 
capacitor  that  changes  its  value  as  a  function  of  both  water  vapor  pressure  and  temperature. 
The  temperature  sensor  uses  a  platinum  100  ohm  resistance  temperature  device  (RTD)  to 
measure  air  temperature. 

The  flux  data  were  stored  in  an  Onset  Model  6  Toshiba  computer  with  500  MB  of 
disk  storage  on  board  the  buoy.  Each  word  is  a  two-byte  data  stored  as  ‘big  endian’ 
hexadecimals  (most  significant  byte  first).  The  mean  meteorology  data  were  stored  in  another 
Onset  Model  6  computer  with  200  MB  of  disk  storage.  The  data  were  stored  in  ASCII 
format.  The  bulk  aerodynamical  method  was  used  to  process  the  mean  meteorological  data. 
C.  ENVIRONMENTAL  INFLUENCE  AND  LIMITATIONS 

The  Monterey  Bay  has  deep  water  with  ocean-like  swells.  The  optical  link  from  Point 
Pinos  to  Marina  State  Beach  crosses  Monterey  Bay  where  there  is  relatively  very  little 
shoreline  influence.  The  buoy  was  located  at  the  mid  point  of  the  link  in  the  deepest  water  and 
•farthest  fi'om  land  'with  open-ocean  ■wind  and  waves. 

The  limitation  in  the  deployment  of  the  buoy  at  a  single  location  leads  to  the 
assumption  of  a  spherically  stratified  atmosphere  "with  a  lateral  homogeneous  refi'activity 
along  the  path  of  propagation.  This  assumption  is  reasonable  only  in  a  well-mixed 
environment  over  the  ocean. 
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In  the  experiment,  the  spatial  average  optical  measurement  is  compared  with  the 
temporal  average  meteorological  data.  In  a  situation  where  there  is  lateral  homogeneity  along 
the  mean  wind  direction,  the  two  averages  could  be  considered  to  be  equivalent  [Ref.  1]. 
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m.  RAY  TRACING 


The  ray  tracing  code  was  used  in  conjunction  with  the  standard  logarithmic 
temperature  profile  and  the  proposed  modified  temperature  profile  to  trace  rays  emanating 
from  the  transmitter  to  the  receiver.  Optical  continuity  is  said  to  be  achieved  if  there  exists 
at  least  one  ray  above  the  ocean  wave  height  connecting  the  transmitter  to  the  receiver. 

The  ray  tracing  technique  used  here  is  based  on  the  methodology  developed  by  Lehn 
[Ref  5],  The  model  assumes  that  the  atmosphere  consists  of  a  number  of  concentric  spherical 
layers  where  the  temperature  is  assumed  to  vary  linearly  within  each  layer.  In  this  case,  the 
curvature  of  a  nearly  horizontal  ray  is  approximately  constant  within  each  layer.  Replacement 
of  circular  arc  by  parabolic  arc  leads  to  a  simpler  quadratic  expression  for  the  ray  trace. 

For  a  ray  traversing  a  stratified  atmosphere  above  a  flat  earth,  Snell’s  Law  dictates 
that 

«sin0  =  A,  (8) 

where  n  is  the  refractive  index  of  the  air  at  any  point  on  the  ray,  0  is  the  angle  between  the 
ray  and  the  vertical,  and  ^  is  a  constant  for  the  ray. 

In  rectangular  coordinates,  the  curvature,  C,  is  given  by 
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(9) 


dh 


C  =  ± 


dx'^ 


1  dz 
1+ — 

3/2 

dx 

where  x  is  the  horizontal  distance  and  z  is  the  elevation  above  the  earth’s  surface.  The  slope 
of  the  ray  is  given  by 

dz  _  cos6 
dx  sin0 

A 


and  the  second  derivative  is  given  by 


dh 

dx^ 


dz  dx 


dz 

dx 


n  dn 


A^  dz 


(11) 


With  Equations  (10)  and  (11),  the  ray  curvature  becomes 


sin0  dn 
n  dz 


(12) 


16 


The  sign  in  Equation  (12)  has  been  chosen  so  that  a  positive  curvature  corresponds  to  a  ray 
that  is  concave  towards  the  earth.  At  optical  frequencies,  the  refractive  index  of  air  can  be 
expressed  as 


n  = 


(13) 


where p  is  the  pressure,  Tis  the  absolute  temperature  in  degrees  Kelvin,  e  =  226x10'®  kg'^ 
and  P  =  3.48x10'^  kg  K.  Note  that  Equation  (13)  is  equivalent  to  Equation  (1)  when  the 
dependence  on  wavelength  is  set  at  A  =  0.83  pm  and  the  pressure  expressed  in  N/m^.  The 
derivative  of  n  is  given  by 


dn 

lb 


€pp 


(14) 


where  the  relation  p  —  pp/Thas  been  used,  p  is  the  density  of  the  air  and  g  is  the  acceleration 
due  to  gravity.  Substituting  Equation  (14)  into  the  expression  for  the  ray  curvature  leads  to 


C  = 

T2 


n 


ep 


P^- 


dz 


(1+ep) 


P^ 


.^L 

dz 


(15) 


17 


In  ray  tracing,  the  circular  arcs  that  represent  the  earth’s  surface,  the  ray  and  the  layer 
boundaries  are  all  approximated  by  parabolic  arcs.  The  coordinate  system  is  shown  in  Figure 
5.  The  surface  of  the  earth  is  given  by 


z  = 


(16) 


where  R  is  the  radius  of  the  earth.  Similarly  a  layer  located  at  a  height  h  above  the  earth’s 
surface  is  given  by 


z  = 


2(R+h) 


+h 


(17) 


An  optical  ray  with  a  radius  curvature,  r  =  1/C,  starting  at  x  =  0,  z  =  and  making  an  angle 
of  (|)o  with  the  horizontal  is  given  by 


x^ 

z  =  -  —  +  xtan(i)(,  +  Z(,  .  (18) 
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Within  the  layer,  the  ray  path  is  given  by  Equation  (18)  with  a  radius  of  curvature 
r,  .  Ray  tracing  is  performed  by  determining  the  intersection  of  the  ray  with  the  layer 
boundaries.  The  7*  layer  lies  betweai  elevations  and  h, .  A  ray  can  enter  a  layer  either  from 
above  or  below  and  could  exit  the  layer  either  from  top  or  bottom  of  the  layer.  Figure  6 
shows  a  ray  entering  the  7'*  layer  from  lower  boundary.  It  will  proceed  to  P  if/*,  >  i?  or  r,  < 
0  and  to  Q  if  r.  is  positive  and  sufficiently  small.  The  coordinates  of  P  are  found  by  solving 
Equations  (17)  and  (18)  with  h  =  A,,  in  Equation  (17)  and  Zg  =  h,.j  in  Equation  (18).  The  x- 
coordinate  of  P  is 


Rr. 


Up 


r-R 


-tan(j)(j  ■ 


<an’<t>o 


-T)\\  1/2 


r,R 


(19) 


If  the  discriminant  is  negative,  the  intersection  P  does  not  exist.  Then  Q  has  to  be  determined. 
This  is  done  by  letting  h  =  /i,.^  in  Equation  (17).  The  x-coordinate  of  Q  is 


Q  ~ 


2i?7-.tan<})Q 


(20) 


Similarly,  for  a  downward-heading  ray,  the  x-coordinate  of  P  at  layer  boundary  hj,j  is  given 


19 


by 


Up  = 


r-R 


-tan(j)o- 


rfi 


1/2 


(21) 


Again,  a  negative  discriminant  implies  that  P  does  not  exist  and  the  ray  bends  upwards  and 
returns  to  layer  boundary  /i,  at  Q.  The  x-coordinate  of  Q  is  determined  from  Equation  (20). 

New  coordinates  are  used  each  time  the  ray  crosses  a  layer  boundary.  Figure  7  shows 
the  case  of  an  upward-heading  ray.  As  the  curvature  is  small,  the  x-coordinate  of  the  new 
origin  is  taken  to  be  Up.  The  layer  boundaries  can  still  be  described  by  Equation  (17).  The 
elevation  angle  of  the  ray  is  given  by 

(j)^'  =  tan‘^/M^  -  tan'^/M^  ,  (22) 


where  and  are  the  ray  and  layer  slopes  respectively  at  P.  The  slopes  are  found  by 
differentiating  Equations  (17)  and  (18)  giving 


tan4)o 


-  m] 

tan(j)o  -  ^  ^  . 


(23) 
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Figure  5.  Ray  Tracing  Coordinate  System  [From  Ref.  5] 


Figure  6.  Upward-heading  Ray  in  Layer  i  [From  Ref  5] 
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Figure  7.  Change  of  Coordinates  when  the  Ray  enters  a  New  Layer  [From  Ref  5] 
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IV.  OCEAN  WAVE  SPECTRAL  ANALYSIS 

A.  MEASUREMENT  APPROACH 

In  order  to  establish  the  possible  effects  of  ocean  surface  waves  on  the  optical  link- 
it  is  necessary  to  examine  the  ocean  surface  wave  spectra  and  wave  heights  at  the  time  when 
each  sample  of  optical  signal  was  measured.  The  wave  information  was  collected  by  the 
motion  pack  and  wave  staff  The  motion  pack  consists  of  three  accelerometers  measuring  the 
accelerations  of  the  x,  y  and  z  directions  (sway,  surge  and  heave)  and  three  rate  gyros 
measuring  the  angular  rates,  0,  ({)  and  t|;  (pitch,  roll  and  yaw).  The  data  was  sampled  at  10.5 
Hz. 

The  true  surface  displacement  is  obtained  by  summing  the  following  three 
components,  viz.,  the  wave  height  as  measured  by  the  wave  staff,  z„,  the  displacement  of  the 
wave  staff  due  to  the  translational  motion  of  the  buoy,  z, ,  and  the  displacement  of  the  wave 
staff  due  to  the  rotational  motion  between  the  wave  staff  and  the  motion  pack  [Ref  6]. 

As  the  wave  staff  and  the  linear  accelerometers  were  mounted  on  the  moving 
reference  of  the  buoy,  it  is  necessary  to  compute  the  surface  displacement  in  an  earth- 
referenced  system.  The  conversion  of  the  buoy-referenced  system  to  the  earth-referenced 
system  can  be  achieved  via  a  transformation  matrix,  Tg^  [Ref  7]  given  by 
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'BE 


COS0  cosijr  sin4)  sin0  cosijr  -  cos<j)  sinijr 
COS0  sini|;  sin(|)  sin0  sim|f  +  cosc})  cosi|r 
-sin0  sin4)cos0 


cos(j)  sin0  cosi|f  +  sincf)  sinijr 
cost})  sin0  sinijr  -  sincf)  cosv|; 

COS(j)COS0 


(24) 


If  ag  represents  the  translational  accelerations  measured  by  the  strap-down 
accelerometers  on  the  buoy,  the  earth-referenced  accelerations,  %  are  given  by 


’^BE  S'  > 


(25) 


where  g  =  (0,  0,  -g)  is  the  acceleration  due  to  gravity.  The  earth-referenced  vertical 
displacement  due  to  the  translational  motion  of  the  buoy  is  given  by 


‘‘lE 


.  (-sin0,cos0sin<f),cos0cos(j>)-g]<iti3& 


(26) 


Similarly,  the  earth-referenced  vertical  displacement  measured  by  the  wave  stalF  is  given  by 


=  ■?;^COS0COS(1)  . 


(27) 
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Finally  the  displacement  due  to  the  rotational  motion  is  determined  by  integrating  the 
rotational  velocity  vector,  Q  x  Tbe  1®  where  Q  is  the  angular  velocity  of  the  buoy  measured 
by  the  rate  gyro  and  Lg  is  the  distance  between  the  accelerometer  and  the  water  surface  at  the 
wave  staff.  The  angular  velocity,  Q  is  given  by 


Q  = 


-6  sim|f  +  <{)  COS0  cosi|r 
6  cosi|r  +  ({)  COS0  sinij; 
t|;-4)sin0 


(28) 


The  vertical  displacement  due  to  the  rotational  motion  is  given  by 


^rE 


-  -  0  sini|;  +  <j)  COS0  cosxjr )  -  0  cost|r  +  ^  COS0  siniff  )]dt  ,  (29) 


where 


’^BE 


(30) 


and  the  subscripts  El  and  E2  represent  the  component  of  the  vector. 
The  true  ocean  surface  displacement  is  given  by 
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(31) 


Once  the  ocean  surface  displacement  is  determined,  the  rms  ocean  surface  displacement,  o 
can  be  computed  and  the  significant  wave  height,  is  given  by 


=  4a  .  (32) 

B.  ANALYSIS  OF  OCEAN  SURFACE  DISPLACEMENT 

The  analysis  of  ocean  surface  displacement  or  elevation  was  carried  in  the  fi'equency 
domain  using  Fourier  Transform.  The  rate  gyros  do  not  have  good  low  frequency  response 
and  therefore  it  is  necessary  to  use  complementary  filtering  to  patch  the  low  frequency  tilt 
reference  of  the  rate  gyros  with  the  linear  accelerometers  as  follows 


0  =  Q^p-sm.-\aylg)^  , 


(33) 


(J)  =  + 


(34) 


and 
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(35) 


ijr  =  ^jjp  +  smr\com)jj,  , 


where  a„  and  represent  the  observed  hnear  acceleration  in  the  x  and  direction,  the 
subscripts  HP  and  LP  represent  the  high  pass  and  low  pass  filtered  output  and  com  is  the 
compass  output.  The  first  tenn  on  the  right  hand  side  of  Equations  (33)  to  (35)  represents  the 
high  pass  angular  motions  fi'om  the  rate  gyros  which  provides  the  wave-induced  angular 
motions  and  the  second  term  represents  the  low  frequency  tilt  reference  fi’om  the 
accelerometers.  The  effect  is  an  all-pass  fiilter  which  removes  the  unwanted  drift  in  the  rate 
gyros  while  retaining  the  low  frequency  tilt  reference. 

To  compute  z  and  ,  the  integration  and  double  integration  are  done  in  the 
frequency  domain  by  dividing  the  Fourier  Transform  of  z  ,.£•  and  z  ^  with  (ico)  and  (ioj)^ 
respectively.  The  division  by  (o  introduces  a  singularity  at  (o  =  0.  However,  there  is  no  ocean 
wave  with  fi-equency  lower  than  0.04  Hz,  i.e.,  25  second  waves,  since  there  is  not  enough 
fetch  in  the  ocean  to  sustain  the  waves  without  breaking.  The  waves  in  the  Atlantic  Ocean  and 
Padfic  Ocean  rarely  exceed  20  seconds  and  25  seconds  respectively  [Ref  8],  The  singularity 
problem  is  alleviated  by  performing  a  high-pass  filtering  of  the  wave  spectrum  at  0.04  Hz.  The 
power  spectral  density  of  the  ocean  surface  waves  and  significant  wave  height  were  computed 
for  the  duration  of  time  when  the  optic  signal  samples  were  being  measured. 
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V.  BULK  AERODYNAMICAL  METHOD 


The  bulk  aerodynamical  formulations  for  the  fluxes  of  momentum,  sensible  heat  and 
latent  heat  across  the  air-sea  interface  and  vertical  profile  estimation  are  given  by 

^  =  PCjjCzXKz)-"/  .  (3(>) 

Qh  =  (37) 

and 

Qe  =  ’  (38) 

where  p  is  the  air  density,  is  the  specific  heat  of  air,  L,  is  the  latent  heat  of  vaporization,  C^, 
Cfj  and  Cg  are  the  height  dependent  bulk  transfer  coefficients  for  momentum,  sensible  heat 
and  latent  heat  respectively,  u,  T,  q  are  the  wind  speed,  potential  temperature  and  specific 
humidity  respectively,  z  is  the  height  above  the  sea  surface  and  the  subscript  s  denotes  the 
respective  values  at  the  sea  surface. 

In  the  atmospheric  surface  these  fluxes  can  be  expressed  according  to  Monin-Obukov 
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Similarity  (MOS)  theory  scaling  parameters  as 


(39) 


Qh  = 


(40) 


and 


Qe  =  > 


(41) 


where  u»,  T*  and  q*  are  the  MOS  surface  layer  scaling  parameters  for  momentum,  temperature 
and  humidity  respectively.  By  combining  Equations  (36)  to  (41),  these  scaling  parameters  can 
be  defined  according  to  the  bulk  transfer  coeflBcients  as 


(42) 


(43) 


and 


(44) 


The  bulk  transfer  coefficients  can  be  defined  as 


Cni?)  = 


(45) 


_ _ 

[ln(2/z J  -  T|;(z/Z)][ln(zV)  -  vlr^r/Z)] 


(46) 


and 


= 


_ _ 

[ln(r/z,„)-tl;(z/Z)][ln(z/zo^)-ilr^(r/i:)]  ’ 


(47) 


where  k  is  the  von  Karman  constant,  and  ijf„,  ilXj.  and  are  the  integrated  forms  of  the 
dimensionless  stability  correction  fimctions  for  wind  speed,  temperature  and  specific  humidity 
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respectively,  given  by 


ilr„(z/I)  = 


2In[(l  +x)/2]  +  ln[(l  +x^)/2]  - 2tan-i(x)  +  ti/2,  z/L<0 
-Sz/L,  z/L^O  ’ 


(48) 


\^t(z/L)  =  ^^(z/L)  = 


21n[(l+x)/2], 

-Iz/L, 


z/L<0 
z/L^O  ’ 


(49) 


X  =  (l-20z/Z)*^'* 


(50) 


and  Zq^  ,  Zqt  and  are  the  roughness  lengths  for  momentum,  temperature  and  humidity 
respectively.  The  Obukov  length  scale,  L  is  defined  as 


L  = 


[K^(7’.+0.6ir^J]  ’ 


(51) 


■where  g  is  the  acceleration  due  to  gra’vity  and  is  the  ■virtual  temperature  introduced  to  take 

into  account  the  effects  of  humidity: 
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r,  =  7(1+0.61^)  . 


(52) 


To  determine  the  vertical  profile  of  wind  speed,  temperature  and  specific  humidity, 
iterative  computation  of  7*,  q*,  Zq^,  Zfyj.,  Zq^  and  L  are  performed  and  the  vertical  profile 
calculated  as 


"  =  w,  +  —  [In  {zizj  - 1|;„  (2/i)]  , 

TC 


(53) 


1C 


(54) 


and 


1C 


(55) 
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The  temperature  profile  estimated  in  this  manner  was  used  as  input  to  the  ray  tracing 
program  to  determine  if  optical  continuity  between  the  transmitter  and  the  receiver  can  be 
estabUshed. 
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VI.  RESULTS  OF  DATA  ANALYSIS 


A.  OCEAN  SURFACE  WAVE  SPECTRA  AND  WAVE  HEIGHTS 

The  power  spectral  densities  of  the  ocean  surface  waves  were  computed  for  the 
duration  when  the  optical  signals  were  being  measured.  A  total  of  41  samples  were  computed. 
A  typical  sample  of  the  wave  spectrum  is  shown  in  Figure  8.  It  can  be  seen  from  the  spectrum 
that  the  ocean  surface  waves  are  dominated  by  swells  with  a  peak  frequency  of  about  0. 1  Hz, 
i.e.  10  second  waves.  There  are  wind  generated  waves  with  a  peak  frequency  at  about  0.4  Hz, 
i.e.,  2.5  second  waves.  The  spectrum  exhibits  a  high  frequency  roll-off  as  o)"^.  This  is  the 
established  behavior  of  the  equilibrium  range  of  wind-sea  spectra  [Ref  9  ]. 

The  rms  wave  height,  o,  was  computed  from  the  power  spectral  density  by  the 
following  relation 


i=l 


(56) 


where  Z(f)  is  the  Fast  Fourier  Transform  (FFT)  of  the  ocean  surface  displacement,  z{t)  and 
Nis  the  sample  size.  For  the  41  samples,  the  rms  wave  height  ranges  from  0. 18  m  to  0.7  m. 
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Figure  8.  Ocean  Surface  Elevation  Spectrum.  The  straight  line  represents  the  -4  slope. 
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B.  TEMPERATURE  SPECTRAL  DENSITY 


The  aspect  of  atmospheric  turbulence  of  relevance  to  optical  propagation  is  the 
variation  of  the  index  of  refraction  or  equivalently,  the  variation  of  the  density  of  air.  Air 
density  depends  on  pressure  and  temperature.  Pressure  fluctuations  are  rapidly  dissipated  and 
their  effects  are  negligible  compared  to  temperature  fluctuations.  Thus,  for  optical 
propagation  purposes,  temperature  fluctuations  are  used  to  characterize  the  optical  effects 
of  turbulence. 

In  the  experiment,  temperature  fluctuations  sampled  at  21  Hz  were  measured  by  the 
sonic  anemometer  at  5.09  meters  above  mean  sea  level.  A  typical  plot  of  the  power  spectral 
density  of  temperature  fluctuations  is  shown  in  Figure  9.  The  spectrum  shows  the  classical 
-5/3  slope  at  high  frequencies.  This  dynamic  condition  in  the  high  frequency  range  agrees 
satirfactorily  with  theoretical  predictions  and  measurements  above  land  [Ref  10].  At  lower 
frequencies,  the  spectrum  shows  evidence  of  ocean  wave  influence.  Figure  10  shows  the 
temperature  spectral  density  and  corresponding  ocean  surface  spectral  density  of  the  wind 
waves.  The  spectral  density  of  wind  waves  was  obtained  by  performing  high  pass  filtering  on 
the  swell-dominated  wave  spectrum,  thus  removing  the  swells  from  the  spectrum.  It  can  be 
seen  that  the  temperature  spectrum  displays  the  characteristic  spectral  density  peak  at  the 
peak  of  the  corresponding  wind  wave  spectrum  indicating  the  influence  of  the  ocean  waves. 
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Figure  9.  Temperature  Spectrum.  The  straight  line  represents  the  -5/3  slope. 
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c. 


OPTICAL  SIGNAL  FLUCTUATIONS 


In  the  analysis  of  the  optical  signals,  the  signal  fluctuations  are  described  by  the  log¬ 
normal  variance  although  the  variance  of  the  signals  can  be  used  as  a  measure  of  fluctuations. 
The  research  does  not  attempt  to  relate  the  log-normal  variance  to  the  saturation 
phenomenon.  The  objective  is  to  analyze  how  the  effects  of  ocean  surface  waves  affect  the 
signal  fluctuations  per  se. 

1.  Log-normal  Distributions 

To  observe  the  signal  fluctuations,  a  histogram  plot  of  the  normalized  log  intensity 
(normalized  to  the  mean  intensity)  was  performed  for  each  signal  sample.  From  the  histogram 
plot,  it  was  found  that  for  signal  samples  where  no  apparent  ocean  wave  occultation  is 
observed,  the  intensity  fluctuations  exhibit  a  log-normal  distribution.  For  signal  samples  where 
wave  occultation  is  evident,  deviation  from  the  log-normal  distribution  results.  Typical 
histogram  plots  of  samples  that  exhibit  log-normal  distribution  and  ‘non  log-normal 
distribution’  are  shown  in  Figures  1 1  and  12  respectively. 

Log-normal  distributed  signals  are  known  to  exist  in  signals  propagating  through 
random  inhomogeneous  medium  [Ref  1 1].  The  distribution  characteristics  are  associated  with 
atmospheric  radio  noise  or  atmospherics  as  it  is  sometimes  called.  It  has  been  found 
empirically  [  Ref  12]  that  the  distribution  P(E>R),  where  E  is  the  instantaneous  value  of  HF 
atmospheric  noise  and  R  is  threshold  voltage  at  the  receiving  point,  approaches  a  Rayleigh 
distribution  for  small  amplitudes.  This  is  due  to  the  interferences  of  many  mutually 
overlapping  small  atmospherics  with  uniformly  distributed  phase.  However  for  large 
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amplitudes,  the  distribution  deviates  considerably  from  the  Rayleigh  distribution  and 
approaches  a  log-normal  distribution.  This  is  believed  to  be  due  to  strong  atmospherics  that 
do  not  overlap  in  time. 

The  theory  of  the  distribution  of  peak  value  of  atmospherics  can  be  explained  as 
follows.  The  peak  value  of  atmospherics,  Ep  can  be  represented  by  the  following  expression 


(57) 


where  Pq  is  the  peak  power  of  the  atmospheric  discharge,  d  is  the  propagation  path,  d.  is  a 
section  of  the  propagation  path,  4  is  the  attenuation  constant  of  the  f  section  and  isT  is  a 
constant.  By  writing  Equation  (57)  as 


=  exp 


(  n 

rS*'*' 


llnPo  -  In^J^.  +  In^ 


(58) 


the  sum  of  the  attenuations  of  each  section.  A,  is  defined.  It  is  assumed  that  the  fluctuations 
oiPo  and  d  are  negligible  compared  to  the  first  term  in  the  exponential  of  Equation  (58). 
Thus,  only  the  first  term  will  thus  be  used  to  determine  the  distribution  of  A.  Since  most 
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atmospherics  are  propagated  from  long  distances  away,  it  is  reasonable  to  assume  that  the 
attenuations  6^  dj  are  independent  and  none  of  the  6^.  is^.’s  make  a  predominant  contribution 
to  the  sum  of  the  attenuations.  Therefore  A  is  normally  distributed  with  probability  density 
function  given  by 


m  = 


exp 


(A-p)^ 


(59) 


2 

where  p  is  the  mean,  is  the  variance.  From  Equation  (58)  we  have 


A  =ln£:^ 


(60) 


and  substituting  Equation  (59)  into  Equation  (60),  the  density  function  o^Ep  is 


KEp)  = 


(Ing^-P)^ 

la\ 


(61) 


Equation  (61)  is  the  log-normal  distribution  which  has  been  found  experimentally  to  apply  to 
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the  propagation  of  optical  signals  through  atmospheric  turbulence  in  the  absence  of 
occultation  of  the  optical  signals  by  ocean  surface  waves. 
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Intensity 


Optical  signal,  30  May  96,  1225  hr 


Minutes 

Optical  signal  distribution 


Figure  11.  Optical  Signal  Intensity  and  Distribution  Plots  and  the  Corresponding  Ocean  Waves, 
in  ocean  waves  plot  indicates  height  of  optical  ray  path  above  the  waves. 
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Optical  signal,  30  May  96,  1859  hr 


Figure  12.  Optical  Intensity  and  Distribution  Plots  and  the  Corresponding  Ocean  Waves  Plot 
(with  Occultation  of  Optical  Signal  by  Ocean  Waves). 
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2. 


Log-normal  Variance  and  Mean 


From  the  above  equations,  the  relationship  between  the  log-normal  mean  intensity  and 
log-normal  variance  can  be  derived.  In  general,  an  electromagnetic  wave  field  E{f,{)  can  be 
represented  by 

E(f,t)  =  exp[7it-f  -  7(0/]  ,  (62) 


where  (o  is  the  temporal  fi’equency  of  the  electromagnetic  field,  k  is  the  propagation  vector 
and  1(f)  is  the  irradiance  or  intensity.  In  the  presence  of  turbulence.  Equation  (62)  becomes 


E(f,t)  =  </(f)>‘^  exp[7^-r  -  7(0/  +  x(^0  + 


where  </(f)>^^^  is  the  average  non  fluctuating  part  of  the  irradiance,  x(r,/)  describes  the 
fluctuating  part  and  <})(r,/)  is  the  phase  variation.  The  random  variables,  x(f,t)  and  ({)(F,/) 
have  Gaussian  probability  distributions.  To  examine  the  irradiance  fluctuations,  the 
instantaneous  irradiance  can  be  written  as 


I(r,t)  =  </(r)>exp[2x(r,/)] 


(64) 


It  follows  that  /(r,  t)  has  a  log-normal  distribution.  Taking  the  average  on  both  sides  of 
Equation  (64),  and  using  the  relation  valid  for  any  real-valued  Gaussian  random  variable. 


{e°^)  =  exp 


a{z)  +  —a^a^ 
2 


(65) 


Equation  (64)  becomes 


</>  =  </>exp(2<x>+2a^^)  . 


(66) 


In  the  presence  of  turbulence,  </>  should  be  a  constant  independent  of  the  turbulence.  It 
follows  that 

<X>  =  .  (67) 


A  plot  of  log  normal  mean  intensity  against  the  log  normal  variance  for  those  signal 
samples  not  occluded  by  the  ocean  surface  waves  in  Figure  13  further  verifies  the  log-normal 
distribution  characteristics  of  the  optical  signals. 
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variance 


Log-nomnal  variance  vs  log-normal  mean 


Log-normal  mean 

Figure  13.  Log-normal  Variance  vs.  Log-normal  Mean. 
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D.  EFFECTS  OF  OCEAN  SURFACE  WAVES  ON  OPTICAL  SIGNAL 
FLUCTUATIONS 

Due  to  the  random  motion  of  the  ocean  waves,  it  is  necessary  to  establish  a  way  to 
measure  the  distance  of  the  ray  path  above  the  ocean  waves.  The  method  adopted  here  is  to 
calculate  the  normalized  mean  square  distance,  h^,  (normalized  to  the  variance  of  the  ocean 
surface  displacement)  of  the  ray  path  above  the  ocean  waves,  given  by 


(68) 


where  o^  is  the  variance  of  the  ocean  surface  displacement,  Zf,,  is  the  lowest  height  of  the 
straight  line  ray  path  above  mean  sea  level,  z,  is  the  instantaneous  ocean  surface  displacement 
and  the  summation  is  performed  over  the  duration  of  the  optical  signal  transmission.  A  large 
mean  square  distance  implies  a  lower  probability  of  the  ray  path  being  blocked  by  the  ocean 
surface  waves.  The  normalized  mean  square  distance  is  plotted  against  the  log-normal 
variance  of  the  signals.  For  sufficiently  large  mean  square  distance,  i.e.,  for  the  mean  square 
distance  greater  than  35,  the  variance  of  the  optical  signal  increases  as  the  mean  square 
distance  decreases.  Those  signal  samples  with  lower  mean  square  distance  also  display  the 
same  trend  but  with  increased  variability,  i.e.,  the  signal  sample  points  are  widely  scattered. 
A  plot  of  log-normal  variance  of  the  optical  signal  against  is  shown  in  Figure  14.  The  signal 
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samples  are  numbered  from  3  to  43  for  easy  identification  of  any  signal  sample.  The  list  of 
optical  signal  samples  is  contained  in  Appendix  A. 

From  Figure  14,  it  can  be  seen  that  there  is  an  apparent  ‘peak’  in  the  signal 
fluctuations  for  signal  samples  14  to  17.  The  higher  signal  fluctuations  are  due  to  the 
increased  roughness  of  the  ocean  wind  waves  associated  with  the  signal  samples.  The  wave 
age,  defined  as  the  ratio  of  the  ocean  surface  waves  phase  speed,  Cp  to  the  wind  speed,  u  is 
a  good  measure  of  roughness  of  the  surface  waves.  Donelan,  et.  al.  [Ref  13],  demonstrated 
that  younger  waves  traveling  slower  than  the  wind  are  rougher  than  mature  waves.  A  plot  of 
the  optical  signal  log-normal  variance  against  wave  age  for  signal  samples  3  to  5  and  14  to 
17  is  shown  in  Figure  15.  The  plot  indicates  a  decrease  in  signal  fluctuations  with  increased 
wave  age.  However,  signal  samples  14  and  15  have  higher  signal  fluctuations  than  samples 
3  and  4  although  they  are  associated  with  wind  generated  waves  of  approximately  the  same 
wave  age.  The  higher  fluctuations  of  samples  14  and  1 5  are  due  to  the  influence  of  the  swell. 
The  influence  of  the  swell  on  sample  14  and  15  is  manifested  in  the  temperature  fluctuations 
spectral  densities  where  the  swell  frequencies  modulate  the  temperature  spectra.  This  effect 
is  shown  in  Figures  16  and  17.  The  corresponding  temperature  spectra  for  samples  3  and  4 
do  not  display  any  prominent  peak  in  the  swell  frequencies.  The  temperature  spectra  and  the 
ocean  surface  elevation  spectra  for  samples  3  and  4  are  shown  in  Figures  18  and  19.  The 
presence  of  swell  caused  a  change  in  the  spectral  distribution  of  the  various  components  of 
the  high  frequency  wind  waves.  There  is  experimental  evidence  that  this  leads  to  an  increased 
in  roughness  of  the  ocean  surface  [Ref  14]. 
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Two  alternative  analyses  were  performed  to  verify  the  trend  of  the  signal  fluctuations. 
In  the  first  method,  it  is  assumed  that  the  ocean  surface  elevation  follows  a  Gaussian 
distribution,  11(0,0).  This  is  true  for  a  sample  with  a  large  number  of  points.  Based  on  the 
height  of  the  ray  path  above  mean  sea  level,  the  probability  of  ocean  surface  waves  protruding 
into  the  ray  path  can  be  computed  fi-om 


Piz>z^)  =  r  -  e  dz  . 


(69) 


The  percentage  of  the  ocean  waves  protruding  into  the  ray  path  is  found  by  multiplying 
Equation  (69)  by  100.  The  percentage  values  are  plotted  against  the  signal  log-normal 
variance.  It  was  found  that  the  signal  fluctuations  exhibit  the  same  trend.  A  plot  of  the 
percentage  of  ocean  waves  protruding  into  the  optical  ray  path  is  shown  in  Figure  20. 

In  the  second  method,  the  log-normal  variance  of  the  signal  samples  are  plotted 
against  the  difierence  between  the  ray  path  height  and  one  significant  wave  height  above  the 
wave  crests.  The  same  trend  in  signal  fluctuations  was  observed.  For  signal  samples  that  are 
less  than  one  significant  wave  height  above  the  ocean  surface  waves,  the  sample  points  are 
widefy  scattered  as  in  the  above  two  analyses.  The  plot  is  shown  in  Figure  21 .  The  analyses 
led  to  the  conclusion  that  the  random  motion  of  the  ocean  waves  is  affecting  air  turbulence 
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immediately  above  the  ocean  surface  and  causing  increased  variability  in  signal  fluctuations 
as  the  ray  path  gets  closer  to  the  ocean  surface. 
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Log-normal  variance 


Log-normal  variance  vs.  h'^2 


Figure  14.  Log-normal  Variance  of  Optical  Signal  vs.  h^. 
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Log-normal  variance 


Log-normal  variance  vs.  Cp/u 


Cp/u 

Figure  15.  Optical  Signal  Log-normal  Variance  vs.  Wind  Wave  Age. 


54 


m''2/Hz  C''2/Hz 


Temperature  spectrum,  31  May  96,  951  hr 


Hz 


Ocean  surface  elevation  soectmm.  31  Mav  96  961  hr 


Figure  16.  Optical  Signal  Sample  No.  14,  Temperature  Spectrum  and  the 
corresponding  Ocean  Surface  Wave  Displacement  Spectrum 
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m''2/Hz  C''2/H2 


Figure  17.  Optical  Signal  Sample  No.  15,  Temperature  Spectrum  and  the 
Corresponding  Ocean  Surface  Wave  Displacement  Spectrum. 
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Ocean  surface  elevation  spectaim,  30  May  96,  1230  hr 

in^  .  .  . . — - - : _ - _ 
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Figure  19.  Optical  Signal  Sample  No.  4,  Temperature  Spectrum  and  the 
corresponding  Ocean  Surface  Wave  Displacement  Spectrum 
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Log-normal  variance 


Log-normal  variance  vs.  percentage  of  ocean  waves  protruding  into  ray  path 


Figure  20.  Percentage  of  Ocean  Waves  Protruding  into  Optical  Ray  Path. 
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Log-normal  variance 


Log-normal  variarKe  vs.  ray  path  1  sig  wave  height  above  ocean  waves 


Figure  21.  Log-normal  Variance  of  Optical  Signal  vs.  Ray  Path  1  Significant  Wave  Height 
above  Ocean  Waves. 


60 


E.  OPTICAL  SIGNAL  MODULATION  BY  OCEAN  SURFACE  WAVES 

The  intensity  spectrum  for  each  optical  signal  sample  was  computed  and  compared 
with  the  corresponding  ocean  surface  wave  spectrum.  The  plot  (Figure  22)  shows  that  for 
optical  signals  that  were  partially  occluded  by  the  ocean  surface  waves,  the  optical  signal 
spectrum  displays  peaks  at  the  swell  frequency  of  the  ocean  waves.  The  optical  signal  is  said 
to  be  modulated  by  the  ocean  waves. 
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Figure  22.  Optical  Signal  Spectrum  and  the  Corresponding  Ocean  Wave  Spectrum. 
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F.  NEAR-SURFACE  MEAN  GRADIENTS 

In  the  above,  only  small-scale  variations  are  considered  whereby  turbulence  produces 
random  fluctuations  in  the  density  of  air  and  hence  fluctuations  in  the  refractive  index.  On 
larger  scale  variations,  the  main  concern  is  on  phenomena  that  affect  the  average  position  of 
the  optical  beam  as  it  traverses  the  atmospheric  medium.  In  the  experiment,  the  effects  of  the 
ocean  waves  on  vertical  temperature  profile  and  hence  the  refractive  index  profile  is 
investigated. 

There  is  evidence  from  the  experiment  that  the  air  turbulence  due  to  the  ocean  waves 
resulted  in  deviations  of  the  vertical  temperature  profile  from  the  logarithmic  profile  thereby 
affecting  the  refraction  of  the  optical  signals.  To  analyze  the  refraction  of  the  optical  signal, 
only  samples  of  those  above  the  ocean  waves,  i.e.,  those  not  occluded  by  ocean  waves,  were 
used.  For  each  time  sample,  the  corresponding  transmitter  and  receiver  heights  above  mean 
sea  levels  were  calculated  from  GPS  surveying  data  and  tide  measurements  provided  by  a 
National  Oceanographic  and  Atmospheric  Administration  tide  gauge.  With  the  measured 
transmitter  and  receiver  heights,  and  using  the  standard  logarithmic  temperature  profile  in  the 
ray  tradng  program,  it  was  observed  that  optical  continuity  between  the  transmitter  and  the 
receiver  could  not  be  established.  An  example  of  the  ray  tracing  plot  is  shown  in  Figure  23. 
The  mechanism  in  which  the  ocean  waves  affect  the  index  of  refraction  needs  to  be  examined. 
This  research  attempts  to  relate  the  effects  of  ocean  waves  on  the  temperature  profile  using 
the  mixing  length  concept. 

The  mixing  length  concept  was  introduced  by  Prandtl  in  1925  to  describe  the 
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properties  of  turbulent  transfer.  In  this  model,  it  is  assumed  that  there  exists  a  length  scale 
called  the  mixing  length  over  which  a  turbulent  eddy  can  travel  and  still  maintain  its  identity. 
Details  of  the  mixing  length  concept  is  contained  in  Appendbc  B. 

The  effects  of  ocean  surfece  waves  on  the  gradients  of  the  vertical  temperature  profile 
are  discussed  in  the  next  section. 
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Height  [m] 


Ray  tracing  using  standard  logarithmic  temperature  profile,  30  May  96,  1225  hr 


Range  [km] 


Figure  23.  Optical  Continuity  cannot  be  established  between  Transmitter  and  Receiver. 
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G,  EFFECTS  OF  OCEAN  SURFACE  WAVES  ON  THE  INDEX  OF 
REFRACTION 

The  refractive  index  at  optical  frequencies  can  be  approximated  by 


n-\  =77.6(1 +7.52x10'^A‘2)(P/7)10‘^  , 


(70) 


where  P  is  the  atmospheric  pressure  in  millibars,  T  is  the  temperature  in  Kelvin  and  A  is  the 
wavelength  of  light  in  pm.  In  the  atmospheric  surface  layer,  pressure  variation  is  small  and 
can  be  neglected.  Variation  in  the  index  of  refraction  is  thus  caused  by  the  variation  of 
temperature  in  the  atmospheric  surface  layer. 

The  vertical  profiles  of  wind  speed,  u,  temperature,  T  and  humidity,  q  all  have 
logarithmic  forms  ^  =(«,  T,  q),  given  by 


In 

/  \ 
z 

'  K 

(71) 


and  the  gradient  profile  given  by 
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(72) 


dz 


where  ^  is  the  flux  scaling  parameter,  is  the  value  at  the  surfece,  z^^^  is  the  roughness  length, 
K  is  the  von  Karman  constant  and  ^^(z/L)is  the  derivative  version  of  the  stability  functions 
given  in  Appendix  B.  In  Equation  (72),  the  term  kz  is  the  turbulence  mixing  length. 

The  above  profiles  are  derived  for  turbulence  above  an  immovable  surface.  A 
distinctive  feature  of  the  boundary  layer  adjacent  to  the  ocean  surface  is  the  presence  of  a 
continuous  random  movement  of  the  ocean  surface.  The  presence  of  the  random  varying  air¬ 
wave  interface  produces  a  constant  variation  in  resistance  of  the  ocean  surface.  The 
roughness  elements  that  originate  from  the  various  components  of  the  wave  spectrum 
contribute  in  varying  degree  to  the  total  wave  dispersion.  The  air  flow  above  the  ocean  waves 
induces  a  pressure  field  which  transfers  part  of  the  energy  of  the  air  flow  to  the  waves  during 
wave  development  or  conversely,  the  waves  transfer  energy  to  the  air  during  swell 
attenuation.  The  surface  waves,  on  the  other  hand,  disturb  the  air  flow  above  them  and  form 
air  waves  which  interact  with  the  gradient  of  the  vertical  profile. 

The  wavy  motion  of  the  surface  aflfects  the  profile  deviation  in  two  ways.  Firstly,  the 
short  waves  (high  fi-equency  component)  contribute  to  the  roughness  of  the  water  surface. 
The  quantity  specifies  the  scale  of  the  turbulent  eddies  that  are  generated  by  the  surface 
roughness  elements.  Experimental  evidence  [Ref  16]  suggested  that  surface  roughness  should 
be  proportional  to  wave  steepness,  H/A,  where  H  is  the  wave  height  and  A  is  the  wavelength 
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of  the  dominant  ocean  wind  waves,  i.e., 


^  ", 
A 


(73) 


As  the  ocean  wave  phase  speed,  C  «  [gA/(2'n:)]^^^,  can  be  expressed  as 


=  a 

Ok 


E. 

2% 


/  \ 


", 

I  C 


(74) 


where  a  is  the  constant  of  proportionality.  Equation  (74)  implies  z^^  «  H  and  z^^  «  w, .  A 
large  u,  implies  a  small  T*  as 


wT 

", 


(75) 


and  hence  a  smaller  temperature  profile  gradient.  In  Equation  (75),  the  overbar  refers  to 
averaging  over  a  suitable  time  period. 

In  this  experiment,  however,  the  wind  speed  was  relatively  low,  with  average  wind 
speed  of  about  3.95  m/s.  The  dependence  of  roughness  on  wave  steepness  may  be  reduced. 
It  may  suflSce  to  use  Chamock’s  formula  that  has  been  interpolated  to  cover  regime  between 
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low  wind  and  high  wind  speed  [Ref.  16], 


=  0.0185 


g 


(76) 


where  g  is  the  gravitational  acceleration  and  u*  is  the  friction  velocity.  The  roughness  length 
for  a  smooth  surface  is  given  by 


(77) 


where  vis  the  dynamic  viscosity  of  air.  The  total  roughness  length  is  given  by 


“Ok 


Z+Z 


(78) 


The  roughness  lengths  z^j.  and  z^^  are  given  by 


^OT  ~ 


0.4- 


(79) 


69 


and 


^0, 


(80) 


Secondly,  the  longer  surface  waves  do  not  contribute  directly  to  the  roughness,  Z(^  but 
they  distort  the  air  flow  immediately  above  the  ocean  surface  and  are  responsible  for  the 
vertical  displacement  of  the  frame  of  reference.  Furthermore,  the  presence  of  the  longer 
waves  causes  a  redistribution  of  wave  energy  and  the  slope  density  of  the  short  waves.  This 
in  turn  influences  the  direction  of  the  mean  stress  vector  by  creating  an  anisotropic  roughness 
field  at  the  surface.  The  long  waves  strained  the  shorter  waves,  causing  them  to  turn  in  the 
direction  of  the  longer  waves.  Donelan  [Ref  17]  showed  that  the  presence  of  long,  gentle 
swell  propagating  in  the  wind  direction  has  a  significant  effect  on  the  wind  sea  spectrum. 
There  is  also  evidence  that  counter  swell,  propagating  against  the  wind  led  to  higher  drag 
coefficients  than  those  ejqperienced  in  pure  wind  waves  [Ref  14],  Therefore  the  presence  of 
swell  will  influence  the  turbulence  of  the  air  above  it.  By  direct  interaction  with  the  air  flow 
and  due  to  the  fact  that  swells  travel  faster  than  the  wind,  the  swells  will  transfer  momentum 
to  the  air  flow  above  it  and  thereby  enhance  turbulent  mixing.  Since  the  random  motion  of  the 
surface  waves  generate  turbulence  in  the  air  above  the  surface  waves  and  increase  turbulent 
mixing,  it  is  appropriate  to  employ  a  statistical  mixing  length  scale  that  accounts  for  the 
increased  mixing.  The  solution  adopted  here  is  to  characterize  the  increased  turbulence  by  an 
effective  random  mixing  length  scale  that  varies  with  the  significant  surface  wave  height, 
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The  significant  wave  height  is  one  of  the  common  parameter  that  oceanographers  use  as  a 
measure  for  the  average  height  of  ocean  surface  waves.  It  is  defined  as  the  average  of  the 
highest  one-third  waves.  The  vertical  profile  of  ^  becomes 


K 


In 


z+H\ 


\  ^0?  ) 


(81) 


and 


dz 


- -  4>£| 


kL 


(82) 


Equation  (82)  indicates  that  the  presence  of  ocean  surface  waves  will  reduce  the  profile 
gradient  in  a  random  manner  that  is  determined  by  the  surface  wave  height.  An  increase  in 
wave  height  will  lead  to  a  reduction  in  the  profile  gradient. 

With  the  modified  temperature  profile,  it  was  demonstrated  by  the  ray  tracing  program 
that  optical  continuity  between  the  transmitter  and  the  receiver  can  be  established.  This  is 
shown  in  Figure  24.  This  approach  is  based  on  applying  the  presumed  scaling  for  the  near 
surface  condition. 
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Ray  tracing  using  modified  temperature  profile,  30  May  96,  1225  hr 


Figure  24.  Ray  Tracing  with  Proposed  Temperature  Profile.  Optical  continuity  is  established. 
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m  CONCLUSIONS  AND  RECOMMENDATIONS 


This  research  showed  that  for  optical  signal  samples  that  are  not  occluded  by  ocean 
surface  waves,  the  signal  fluctuations  exhibit  a  log-normal  distribution.  This  is  characteristic 
of  signals  propagating  through  a  random  inhomogeneous  medium.  For  signal  samples  that  are 
occasionally  obscured  by  ocean  surface  waves,  the  signal  fluctuations  deviate  fi-om  a  log¬ 
normal  distribution. 

The  effects  of  ocean  surface  waves  on  optical  signal  fluctuations  are  revealed  in  the 
temperature  fluctuation  spectral  densities.  The  temperature  spectral  densities  display 
prominent  peaks  at  the  corresponding  peak  frequencies  of  the  ocean  surface  waves.  These 
characteristic  peaks  are  not  present  in  the  temperature  spectral  densities  computed  fi'om 
meteorological  measurements  over  land. 

Optical  signal  fluctuations  increase  as  the  signal  ray  path  gets  closer  to  the  ocean 
surface  waves  demonstrating  the  sigmficant  effects  of  the  ocean  surface  waves  on  signal 
fluctuations.  Signal  fluctuations  were  also  shown  to  be  dependent  on  wave  age  and  the 
influence  of  ocean  swells.  Younger  waves  are  rougher  and  are  likely  to  cause  a  greater 
increase  in  agnal  fluctuations  than  more  mature  waves  cause.  When  heights  of  signal  samples 
are  sufficiently  low  so  that  the  ocean  waves  occasionally  obstruct  the  ray  paths,  the  ocean 
waves  modulate  the  optical  signals.  Under  moderate  or  high  wind  conditions,  the  ocean  wind 
wave  steepness  should  be  incorporated  into  the  roughness  length  scale,  . 

On  larger  scale  variations,  the  ocean  surface  waves  modify  the  vertical  temperature 


73 


profile  and  hence  the  refi’active  index  profile.  The  standard  logarithmic  temperature  profile 
derived  for  meteorological  conditions  over  land  is  not  suitable  for  the  marine  environment 
above  the  undulating  ocean  surface.  Ray  tracing  using  the  standard  logarithmic  temperature 
profile  indicates  that  optical  continuity  could  not  be  established  between  the  transmitter  and 
the  receiver.  Unlike  conditions  over  land,  the  ocean  waves  distort  the  air  flow  above  and 
increase  turbulent  mixing.  This  increased  mixing  is  proportional  to  the  ocean  wave  height. 
The  increased  turbulent  mixing  can  be  accounted  for  by  incorporating  a  statistical  mixing 
length  proportional  to  the  ocean  wave  height.  The  proposed  temperature  profile  assures 
optical  continuity  between  the  transmitter  and  the  receiver. 

The  research  has  demonstrated  the  significant  effects  of  ocean  surface  waves  on  the 
propagation  of  low  altitude  optical  signals  above  the  ocean.  The  results  of  this  study  provide 
an  insight  into  the  effects  of  ocean  surface  waves  on  an  optical  signal  skimming  the  ocean 
surface  from  an  engineering  view  point.  The  study  is  not  intended  to  disprove  existing 
theories  and  hypotheses  but  rather  to  emphasize  the  relative  importance  of  certain  parameters 
that  affect  the  propagation  of  optical  signals  in  the  vicinity  of  the  ocean  surface  where  the 
influence  of  ocean  surface  waves  is  significant. 

The  easting  model  (e.g.  IR  Tool)  used  by  the  Navy  utilizes  the  standard  logarithmic 
temperature  profile  to  predict  refraction  of  signals.  The  standard  logarithmic  profile  is  only 
valid  for  heights  that  are  several  significant  wave  heights  above  the  ocean  surface  waves 
where  the  ocean  wave  effects  are  negligible.  Within  a  few  meters  of  the  ocean  surface  waves, 
excessive  ray  bending  occurs  due  to  the  logarithmic  singularity  of  the  temperature  profile  as 
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the  ray  approaches  the  ocean  surface.  The  proposed  statistical  mixing  length  approach  to 
modifying  the  temperature  profile  is  perceived  as  an  engineering  fix  to  alleviate  the  large 
gradient  in  this  critical  layer  above  the  ocean  surface  where  sea-skimming  threats  are 
predominant.  Further  refinement  in  the  modeling  is  recommended  through  extensive 
experimental  measurements  and  data  analysis.  Modeling  of  signal  propagation  above  an  air- 
ocean  interface  spans  an  inter-disdplinary  study  of  electronic  engineering,  oceanography  and 
meteorology.  The  complexities  of  the  air-ocean  interface  presents  a  unique  challenge  to 
scientists  and  engineers  in  a  continuous  search  for  more  accurate  modeling  of  signal 
propagation  in  this  critical  marine  layer. 
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APPENDIX  A.  OPTICAL  SIGNAL  SAMPLES 
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Optical  Signal  Samples 


B  C  D  E 


Date  time  Rx  ht  [m]  Txht[m]  height  of  ims  % 


_ ray  path  wave  ht _ 

5301225  4.69  5.16  1.71  0.36  0.000106 

5301230  4.56  5.18  1.64  0.36  0.000248 


1.789256 


18.7128 

16.5491 


84 

_ 68 

10.6063 1 


H 


log  var 


0.079 

0.082 


95 

67 


log  var 


[wave  occ) 


843 

4.26 

859 

3.86 

.36 

0.58 

0.956063 

7.9751 

.11 

0.45 

0.688336 

5.7584 

1.123754 


5.1516 


.8771 

.6558 


0.19 

3.23E-11 

39.5961 

0.065 

1.22E-07 

1.06E-05 


36.3271 

40.2262 


0.123 

0.108 


I  1.27 

0.23 

1.88E-06 

37.5196 

0.2 

4.1E-08 

35.2933 

0.33 

0.030874 

15.0818 

5311834 

3.69 

5311836 

3.53 

5311839 

3.43 

5311932 

2.94 

5311937 

3.07 

5311948 

3.32 

6012208 

3.22 

0.06726 

0.016937 

0.50895 


11.8792 

9.2211 

8.2193 


ESil 


6012243 

6012308 

6021245 


6021417 

6021451 


6030930 


6030937 

6031540 

6031551 

6031605 

6031611 

6041325 

6041417 

6041426 

6051933 

6051945 

6052000 

6052041 


3.5517 


1  0.2 

2.511705 

3.3162 

0.252 

’  0.19 

0.134146 

6.4756 

0.134 

0.18 

0.010306 

12.4613 

0.257 

■■■ 

85.2618 
99.52231  " 
0.015119’ 
0.117065’ 
0.004213' 
0.001032’ 
0.052834  ’ 
0.079321 ~ 
0.020864 
9.854486  ’ 
0.249417’ 
0.703241 
20.57304 
99.9921 
99.99999 


3.0506 
7.2935 
9.949 
6.8136 
24.3964 
17.7704 
11.5735 
13.1283 
11.8525 
3.1197 
7.3583 ' 
10.8039" 
1.5127 
8.4996 ' 
19.4317" 


0.107 

0.126 

0.103 

0.091 


0.157 

■0?i45 

0.181 


0.132 

0.079 

0.269 


0.196 

0.123 

0.077 
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Legend: 
Column  A; 
Colunm  B: 
Column  C: 
Column  D; 
Column  E; 
Column  F: 
Column  G: 

Column  H: 
Column  I; 


Date  and  time  of  signal  sample. 

Receiver  height  above  mean  sea  level  (in  meters). 

Transmitter  height  above  mean  sea  level  (in  meters). 

Lowest  height  of  the  ray  path  above  mean  sea  level  (in  meters). 

Rms  ocean  wave  height  (in  meters). 

Percentage  of  ocean  surface  waves  protruding  into  optical  signal  ray  path. 
Normalized  mean  square  distance  of  optical  signal  ray  path  above  ocean 
surface  waves. 

Log-normal  variance  of  signal  samples  not  occluded  by  ocean  surface  waves. 
Log-normal  variance  of  signal  samples  occluded  by  ocean  surface  waves.  This 
is  not  strictly  correct  as  the  signal  samples  no  longer  exhibit  log-normal 
distribution. 
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APPENDIX  B.  MIXING  LENGTH  HYPOTHESIS 

The  mixing  length  concept  was  introduced  by  Prandtl  in  1925  to  describe  the 
properties  of  turbulent  transfer.  In  this  model,  it  is  assumed  that  there  exists  a  length  scale 
called  the  mbdng  length  over  which  a  turbulent  eddy  can  travel  and  still  maintain  its  identity. 

Suppose,  under  neutral  atmospheric  conditions,  a  parcel  of  air  moves  vertically 
through  a  small  distance,  /,  from  its  original  position,  z+l,  where  the  characteristic  quantity, 
q(z),  has  an  average  value  of  q{z+T),  then 

q{z)  =  q(z+I)  .  (83) 

Using  R^old’s  decomposition,  q(z)  can  be  written  as  a  sum  of  an  average  component  and 
a  fluctuating  component, 

qiz)  =  q(z)  +  q'(z)  . 

Using  mixing  length  hypothesis, 

q'(z)  =  q{z+t)  -  q{z)  »  M  . 

dz 


(84) 


(85) 


Suppose  —  >  0,  then  as  shown  in  Figure  25,  /  >  0  implies  that  the  eddy  comes  from 
dz 

above  and  hence  the  vertical  turbulent  velocity,  w'  <  0 .  By  the  same  token,  /  <  0  implies 
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w'  >  0.  The  above  argument  indicates  that  / and  w'  are  negatively  correlated,  i.e.. 


w'l  <  0  . 


(86) 


Assuming  that  mixing  length  hypothesis  applies  to  the  turbulent  velocity  component 
and  the  air  temperature  then 


TXz)  =  I— 
dz 


(87) 


and 


»'(z)  =  . 

dz 


(88) 


where  T'  is  the  turbulent  air  temperature  and  u'  is  the  turbulent  streamline  wind  velocity. 
Equation  (88)  implies 


w’u'  =  w'l  — 
dz 


(89) 


du 

In  the  atmospheric  surface  layer,  —  >0  and  since  w'l  <  0,  then 

dz 
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w'u'  <  0  . 


(90) 


Assuming  the  eddy  velocity  components,  m'  and  have  the  same  order  of  magnitude,  then. 


w 


u 


lis. 

dz 


(91) 


Equations  (89),  (90)  and  (91)  lead  to 


(92) 


In  the  atmospheric  surface  layer,  the  mixing  length  varies  with  height  above  the  earth’s 
surface.  Near  the  earth’s  surface,  mixing  cannot  take  place  over  large  distances  due  to  the 
solid  boundary.  Away  from  the  surface,  the  turbulence  is  less  inhibited  and  miving  is 
enhanced.  In  general,  it  is  assumed  that 

/  =  Kz  ,  (93) 


where  k  is  the  von  Karman  constant. 

In  the  atmospheric  surface  layer,  the  main  task  is  to  establish  relationships  between 
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vertical  gradient  of  mean  quantities  of  wind  speed,  temperature  and  specific  humidity  and 
surface  fluxes.  In  particular,  the  surface  fluxes  of  momentum,  t,  sensible  heat,  Q^,  and  latent 
heat,  Qe,  are  given  by 


T  = 


T~7  2 

-puw  =  , 


(94) 


a 


H 


=  pc^Tw'  = 


(95) 


and 


Qe  =  pL^q'w'  =  -pL^u^q^  . 


(96) 


where  is  the  specific  heat  of  air,  and  is  the  latent  heat  of  vaporization  of  water. 
Combining  Equations  (92),  (93)  and  (94), 


^  =  ill 

dz  KZ 


(97) 


where  u  is  the  mean  wind  velocity  (the  overbar  has  been  dropped  for  convenience).  Similar 
equations  for  T  and  q  can  be  written; 
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^  =  ii 

dz  KZ 


(98) 


and 


dz  KZ 


(99) 


For  non-neutral  atmospheric  conditions.  Equations  (97),  (98)  and  (99)  become 


du 

dz 


KZ 


(100) 


dz 


=  ^  <j,^z/L) 

KZ 


(101) 


and 
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where  (j)’s  are  the  stability  functions  determined  empirically, 


U 


(1  zlL<0 

\+%zlL,  zlL>0  ’ 


(103) 


4>r  =  4)^ 


(1  -  \6zlLy^'^, 

\+lzlL, 


zlL<0 

zlL>0 


(104) 


and  z/L  <  0,  z/L  >0,  and  z/L  =  0  correspond  to  unstable,  stable  and  neutral  stratification 
respectively  and  L  is  the  Obukhov  length. 

Integrating  equations  (100),  (101)  and  (102)  yield  profiles  for  mean  u,  T  and  q, 


u=u^^  —  [\n(zlzJ-}if^{zlD)\  , 

K 


(105) 


(106) 


T=T^^. 


[ln(z/2„„)-t|f„(z/L)] 


and 


+  — [ln(z/Zop-i|;^(z/Z)]  , 
K 


(107) 


where  ,  T, ,  q  axt  the  values  at  the  surface  and  i|;’s  are  the  empirically  derived  stability 
correction  functions  given  in  Chapter  V. 
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Figure  25.  Mixing  Length 
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